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ABSTRACT 

Despite their importance, very few RR Lyrae (RRL) stars have been known to reside 
in binary systems. We report on a search for binary RRL in the OGLE-III Galactic 
bulge data. Our approach consists in the search for evidence of the light-travel time 
effect in so-called observed minus calculated (O — C) diagrams. Analysis of 1952 well- 
observed fundamental-mode RRL in the OGLE-III data revealed an initial sample of 
29 candidates. We used the recently released OGLE-IV data to extend the baselines up 
to 17 years, leading to a final sample of 12 firm binary candidates. We provide O — C 
diagrams and binary parameters for this final sample, and also discuss the properties 
of 8 additional candidate binaries whose parameters cannot be firmly determined at 
present. We also estimate that > 4 per cent of the RRL reside in binary systems. 

Key words: binaries: general - methods: data analysis - stars: oscillations - stars: 
variables: RR Lyrae - techniques: photometric - stars: fundamental parameters. 


1 INTRODUCTION 


RR Lyrae (RRL) stars play a key role in astrophysics. They 
are important distance indicators, allowing us to determine 
the distance to th e closest galaxies (e.g., ICacciar 3 120131: 
iDambis et al.l [20131 ) . and thus providing an important step 
in the calibration of the extragalactic distance scale. Their 
importance in the context of galaxy form ation and evolu tion 
is also being increasingly recognized (e.g.. lCatelanll2009l) . In¬ 
deed, RRL stars being unmistakably old, their distribution 
in the Galactic halo pro vides evidence of th e early Milky 
Way formation history llDrake et al.l l2013al lbl: ISesar et al.l 
120131 : iTorrealba et alJ l2015h . In addition. RRL stars help 
trace the spatial distribution, and even the ag e, of some 
of the Milky Way’s oldest stellar populations (Eeeli 19921; 
ICatelan fe de Freitas Pachecdl 19931 : IP Aanv et al. 2013l b 

Pespite their usefulness, there are still important un¬ 
certainties affecting the fundamental parameters and phys¬ 
ical properties of these stars. Trigonometric parallaxes of 
even the closest RRL are notorious for their large error bars 
dBenedict et al.ll201ll ). To complicate matters further, the 
exact evolutionary status of even RR Lyr itself is uncertain, 
with indications that it may be significantly overluminous, 
compared to th e zero-age horizontal branch (ZAHB) level at 
its metallicity dCatelan fe Cortisll2008l : iFeast et al.ll2008l ). 

RRL variables occupy a fairly narrow strip - the so- 
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called instability strip - at intermediate temperatures along 
the HB. As such, their exact mass value is crucial in estab¬ 
lishing whether an HB star will ever become an RRL, or 
instead, a non-variable blue or red HB star. In this sense, 
theory predicts that the masses of RRL stars should decrease 
with increas ing metallicity , with little sc atter at any given 
[Fe/H] (e.g.. ICatelanlll992 ; ISandagdl2006l ). Pirect empirical 
confirmation of this important result is, however, still lack¬ 
ing. Most of the available information regardi ng RRL masses 
comes from the so-called Petersen diagram (IPetersenll 19731 ) 
of double-mode RRL (RRd) stars. RRd stars are observed 
to pulsate simultaneously in the fundamental and first over¬ 
tone radial modes, and their distribution in the period ratio 
vs. period (Petersen) diagram is predicted to be a strong 
function of the pulsating star’s mas s, in addition to othe r 
parameters, such as metallicity (e.g.. IPopielski et al.ll200(il ). 
RRL star masses can thus be derived by comparing the ob¬ 
served positions of RRd stars in the Petersen diagram with 
those predic ted according to stellar evolution and pulsation 
theory (e.g., iBono et~aDl 19961 : IPekanv et al.l2008l ). However, 
this method can only be trusted to provide accurate masses 
if the theoretical framework upon which it is based is itself 
accurate. To constrain the theories themselves, it is imper¬ 
ative to obtain a model-independent mass measurement. 

Binary systems allow the derivation of the masses of 
its components, if the orbital parameters are known. In 
the case of binary systems containing classical Cepheids, 
analysis of their orbital parameters has played a crucial 
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role in accurately esta blishing their physical properties, 
including their masses (IPietrzvriski et all 120101 1 . Further¬ 
more, Cejaheid-bearing binary systems are relatively com¬ 
mon (lSzabadoal2003l l ■ and even systems in wh ich both com¬ 
pone nts are Cepheids are now known to exist (iGieren et al.l 
1201411 . 

The situation regarding RRL stars could hardly be more 
different, as for long only one RRL, TU UMa, has been 
kno wn to reside in a b inary system (see lSaha fe Whit3ll990l 
and lWade et al.lll999l . for very detailed analyses of this star). 
Recently, eclipsing binary RRL candidates have been found 
by th e OGLE project in the G alactic bulge IS oszyns ki et al.l 
EB. as well as in the LMC (ISoszvriski et al.ll2009l h How- 
ever, follow-up observations and modeling of the bulge can¬ 
didate have shown that the pulsating component of this 
system has too low a mass to be a bona -fide RRL star 
dPietrzvriski et akl 120121 : ISmolec et al.l l2013l l . Other candi¬ 
date e clipsing RRL found b y the OGLE-II project in the 
LMC ((Sosz^ns^et^lJbOOjj) have turned out to be optical 
blends dPrsa et~ ImT Careful analyses of the RRL LCs 
of the Kepler mi ssion have uncovered only three possible bi¬ 
nary candidates dLi fe Qianll20l3 : lGuggenberger fe Steixneil 

B) . 

Recently, a number of metal-poor, car bon-rich RRLs 
have been identified via spectroscopy (e.g., iKennedv et al.l 
120141 1. The anomalous C abundance of these stars can be 
explained by mass accretion from a more massive compan- 
ion, which has evolved through the asymptotic giant phase 
dStancliffe et al.ll2013l l. suggesting that these RRLs might 
be members of binary systems in which the other compo¬ 
nent has already evolved to the white dwarf stage. However, 
the binary nature of these RRLs has not yet been directly 
established. 

I 11 this Letter , we describe our search for pe riodic 
phase variations caused by the light-travel time effect (llrwinl 
EH) in the LCs of a subsample of well-observed RRL 
stars from the OGLE- III survey towards the Galactic bulge 
dSoszvriski et al.1 1201 lh . We augment the LCs of the binary 
candidates thus obtained wi th newly published ph otometry 
from the OGLE-IV survey (ISoszynski et al.ll2014h . in order 
to increase the observational baseline. We derive binary pa¬ 
rameters for the best candidates. Furthermore, we discuss 
the RRL binary fraction, as well as the detectability of such 
systems through the light-travel time effect. 


2 DATA AND ANALYSIS 

Our initial an alysis is based on the /-band OGLE-III LCs 
for bulge RRL (ISoszynski et al.ll201ll h We analyze stars with 
fundamental-mode pulsation (RRab subtype), as well as 
having an observational baseline covering more than 10 yr. 
These criteria define a subsample of 1952 RRab variables. 
For these, we have utilized the s o-called observ ed minus cal¬ 
culated (O — C ) diagrams (e.g., ISterkenll2005l l. adopting a 
linear ephemeris (C) for the variables: 

C(t) — to + /puis E, (1) 

where to is the initial epoch, P pu \ s is the pulsation period, 
and E is the epoch number, corresponding to the number 
of elapsed pulsation cycles since to- Times have been trans¬ 
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Figure 1 . O — C analysis of the OGLE-III observations, (a) 
The original LC folded with the pulsation period ( grey dots), 
its Fourier fit ( red line), one section of the LC ( black dots), and 
the least-squares template fit ( blue dashed line), (b) The result¬ 
ing O — C diagram (dots), and the corresponding binary-model fit 
(blue line). The point corresponding to the measurement shown 
in (b) is indicated by an arrow, (c) The folded LC, after correct¬ 
ing the times of observations for the binary motion (panel [b]). 
The Fourier fit ( blue dashed line) now follows the real LC shape 
much more closely than the original fit (red line). 


formed into Barycentric Julian Dates (BJD) in the Ba rycen- 
tric Dynamical Time standard (lEastman et al.|[20ld l. 

In the case of variable stars, the O — C diagram is most 
commonly constructed by subtracting the ephemeris (C) 
from timing observations ( O ) of particular features of the 
LCs, such as maxima or minima, a nd plotting this quantity 
as a function of time. [Hertzsprunel (jl919l l proposed utilizing 
the whole LC for deriving O — C points of a variable by fit¬ 
ting the phase of a LC template to the observations. Due to 
the generally sparse sampling of OGLE observations, deter¬ 
mining the times of indivi dual ma xima at different epochs 
is unfeasible, making iHertzsprunel ’s method vastly superior 
for obtaining the required O — C measurements. We thus 
adopt the latter method in our analysis. 
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Figure |T] illustrates the procedure for one of the binary 
candidates. First we created a LC template b y fitting th e 
original LC with a Fourier series using lcfit (ISodoill2012lh 
as shown in Figure [ljt. The LCs were divided into short sec¬ 
tions corresponding to different observing seasons. Sections 
longer than 160 days were split in two, in order to achieve 
better time resolution. We have derived the O — C points 
by least-squares fitting the LC templates in phase to each of 
these segments. These phase shifts were then used to con¬ 
struct the O — C diagrams (Fig. |T}z>). 

2.1 Selection of binary candidates 

The light-travel time effect manifests itself as a strictly peri¬ 
odic phase modulation of the LC of a variable star, caused by 
its orbital motion around the common center of mass in a bi¬ 
nary system. The change in the observed times of particular 
features of the LC, and consequently in the corresponding 
O — C values, have the form: 



1 + ecos(i') 


where a is the semi-major axis, i is the inclination, e is the 
eccentricity, and to is the argument of the periastron. The 
true anomaly v is a function of the time t, the orbital period 
-P or b, the time of periastron passage T per i, and e. 

We have visually inspected the O — C diagrams con¬ 
structed as described previously, in search of O — C shapes 
allowed by Eg. l2l fe.g.. Ilrwinlll959l b In some cases, particu - 
larly when the period values found bv ISoszvriski et all (1201 il l 
have turned out to be inadequate, we have constructed the 
O — C diagrams with several different pulsation periods. 

Figure [0 shows one of the best candidates for a binary 
system containing an RRL star obtained in this way. In panel 
b of this figure, a preliminary orbital solution is also pro¬ 
vided. Figure |TJ: demonstrates the accuracy of the fit; the 
prewhitening for the orbital solution results in a LC which 
has an rms deviation of only 12 mmag. 

A large fraction of stars show errati c period chan ges, 
poss ibly connected to t he Blazhko effect jBIaz kcil 119071 : see 
also ljurcsik et al.[l201ll l. The Blazhko effect manifests itself 
as an amplitude and phase modulation of the LC of RRL 
stars, and accordingly, variables affected by it have unreli¬ 
able O — C diagrams for our purposes. For this reason, we 
have discarded binary candidates for which the LC shape 
and/or amplitude showed correlated changes with the phase 
variations, which might be indicative of the Blazhko effect. 
However, the presence of the Blazhko effect is unclear in 
some of the variables, as blending with another variable star 
and/or large photometric errors (in case of the fainter stars 
in the subsample), together with the sparse coverage of the 
OGLE-III data, could mimic and/or mask the change of the 
LC shape/amplitude of Blazhko variables. 

2.2 Analysis of the candidates 

Following the above procedure, we have selected 29 poten¬ 
tial binary candidates. In order to clarify the status of the 
stars in this sample, and also to refine the binary parameters 
of the candidates, we next combined the OGLE-III and IV 
datasets for these stars, thus increasing the baselines of their 


Table 1. OGLE IDs“ of Uncertain RRL Binary Candidates 


Likely non-binaries 

Likely binaries with poor fits 

Symptom 

ID 

Symptom 

ID 

Blazhko effect 

05075 

P or b > 6000 d 

04522 


07773 


05135 


10519 


05152 


12027 


10891 


13260 


11683 


13698 


12611 




12698 

Period changes 

05778 

Very small amplitude 

14852 


06876 




12195 




“OGLE-BLG-RRLYR 


O — C diagrams. The full analysis of the combined OGLE-III 
and IV datasets will be reported in a future paper. 

We have combined the OGLE-III and IV datasets for 
each of our 29 candidates by correcting for t he difference be¬ 
tween their average magnitu d es, as given bv ISoszvriski et all 
d2009lf and ISoszvriski et alJ d20l4l . We then repeated the 
previous analysis on the data thus combined. We have in¬ 
spected the new O — C diagrams and the combined LCs, 
and have found strong evidence that 6 of the candidates do 
in fact show the Blazhko effect, which was not evident from 
the relatively sparse OGLE-III data alone. The O — C varia¬ 
tions of two additional candidates appear to be caused com¬ 
pletely by their linear period changes, while for a third star, 
it is probably caused by irregular period changes. These dis¬ 
carded candidates are listed in the second column of TableQ] 
leaving us with a refined sample of 20 binary candidates. 

Eq. [5] does not take possible changes in the pulsation 
periods of the RRL stars into account. Such changes can 
significantly alter the O — C diagram, even on relatively 
short (~ 10 yr) timescales. If the period changes linearly 
with time, the O — C diagram has the following shape: 

(O — C){t) = co + cR + C2i 2 , (3) 

where C 2 = /3 is the linear period-change rate. The increased 
time-base of the combined datasets allows us to take this 
effect fully into account. We thus fit the O — C diagrams of 
the 20 remaining candidates with the sum of Eqs. [2] and [3] 
and find that reliable parameters can be derived for 12 of 
them. The other 8 variables either have orbital periods that 
are presumably comparable or longer than the baseline of 
the available observations (therefore the parameters of the 
fit are degenerate), or have very small O — C amplitudes. 
These stars, which may still be bona-fide binaries and thus 
merit continued monitoring, are also listed in Table Q] 

In order to derive the best possible parameters for this 
12-star sample, we iterate the O — C solution once, by means 
of an improved Fourier fit to a (O — C)-subtracted LC. We 
also 3cr-clip the LCs, discarding the worst-quality (typically 
~ 1 per cent) observations for each star. The final O — C 
points are determined by fitting these templates to the LC 
segments. The final diagrams, together with their fits, are 
shown in Figure [2] Table [2]gives the fitted binary parameters 
for each star, as well as other relevant parameters of the fit. 


© 2015 RAS, MNRAS 000, □[6] 





























4 G. Hajdu et al. 


07640 P pll | s =0.5537702 d P orP =1255d 14408 P pu | s =0.4460234 d P prP =2734 d 04376 P pu , s =0.4910175 d P orb =27S2 d 


0.00 

g 

o 

6 

-o.oi 


- 0.02 



0.00 

g 

o 

6 

- 0.01 





0 1000 2000 3000 4000 5000 6000 7000 0 1000 2000 3000 4000 5000 6000 7000 0 1000 2000 3000 4000 5000 6000 7000 

BJD-2450000 BJD-2450000 BJD-2450000 


Figure 2 .O — C diagrams of RRL binary candidates. The O—C points (filled circles-. OGLE-III data; empty circles: OGLE-IV data) 
are fitted with the sum of a linear period change (Eq.[3] dashed line) and the binary orbit (Eq.(2] dot-dashed line, shifted downwards for 
clarity). For each star, the OGLE ID, pulsation period, and orbital period are given on top of the upper x-axis. 


3 DISCUSSION 

We have completed the first systematic search for binaries 
in a subsample of OGLE Galactic bulge RRL stars utilizing 
the light-travel time effect. Twenty probable binaries have 
been found analyzing the O — C diagrams and LCs of 1952 
OGLE-III bulge RRab variables, which represents about 1 
per cent of this particular subsample (fundamental-mode 
pulsation, >10 yr observational baseline). This allows us to 
very roughly constrain the RRL binary fraction, as follows. 

Approxima tely 50 per cent o f RRab stars show th e 
Blazhko effect djurcsik et al.l 120091 ; iKolenberg et al.1 l2010h . 
Due to their LC phase and shape changes, determining their 
binarity through the O — C method is impossible. Therefore, 
all of the Blazhko binaries have necessarily been missed, and 
so the total binary fraction must be closer to 2 per cent than 
to 1 per cent. We assume that the fraction of stars whose bi¬ 
narity was missed due to erratic pulsation period changes is 


small in comparison. Binaries with very low inclinations and 
very long periods are also missed. As we have no informa¬ 
tion about the fraction of long-period binaries in the sample, 
and as our method is more sensitive to stars with high in¬ 
clination (because the O — C amplitude is proportional to 
asini), we conservatively assume that at least half of the 
binaries are still missed due to these two selection effects. 
Based on these arguments, presumably at least 4 per cent of 
RRL variables reside in binary systems in the sample. The 
recovery rate using the O — C method is thus ~ 25 per cent, 
or perhaps lower. 

Very few RR. Lyrae stars have been reported to reside 
in binary systems yet (see Sect. [I]). The current sample al¬ 
lows us to assess the observational requirements to discover 
additional such systems. The period distribution of the bi¬ 
nary candidates is highly skewed: we find no binaries with 
orbital periods shorter than ~ 3.5 yr, and there is a clus- 
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Table 2. Fitted and derived parameters of the RRL binary candidates 11 


ID 

-^orb 

(d) 

T ■ 

peri 

(d) 

e 

L J 

(deg) 

a sin i 

(AU) 

a 

(AU) 

P 

(d Myr" 1 ) 

K 

(km s _1 ) 

/M 

(M 0 ) 

07640 

1255 ± 4 

9598 ± 54 

0.16 ±0.04 

-30 ± 15 

0.97 ± 0.02 

0.16 

0.05 ±0.01 

8.5 

0.0774 

14408 

2734 ± 38 

8871 ± 78 

0.60 ±0.12 

162 ±8 

0.30 ±0.03 

0.11 

0.16 ±0.01 

1.5 

0.0005 

04376 

2782 ± 19 

8147 ± 53 

0.79 ±0.11 

-92 ±6 

0.74 ± 0.06 

0.11 

0.04 ±0.01 

4.7 

0.0070 

06498 

2789 ±18 

8137 ± 134 

0.12 ±0.04 

-82 ± 16 

2.35 ±0.05 

0.16 

-0.11 ± 0.03 

9.2 

0.2228 

07566 

3485 ± 15 

8612 ±33 

0.54 ±0.03 

-178 ±3 

1.72 ±0.04 

0.14 

-0.03 ±0.02 

6.4 

0.0556 

09789 

3674 ± 30 

7221 ±81 

0.18 ±0.04 

78 ± 10 

2.80 ±0.06 

0.23 

0.24 ±0.03 

8.4 

0.2162 

13534 

3715 ± 27 

6327 ± 73 

0.26 ±0.03 

20 ± 5 

1.00 ±0.02 

0.22 

0.02 ±0.01 

3.0 

0.0096 

07943 

3726 ± 40 

9273 ± 150 

0.14 ±0.03 

-15 ± 13 

1.28 ±0.03 

0.20 

-0.05 ±0.02 

3.8 

0.0200 

11522 

3894 ± 153 

8751± 411 

0.30 ±0.14 

-61 ± 24 

0.57 ± 0.08 

0.12 

0.37 ± 0.04 

1.7 

0.0017 

06981 

4014 ± 85 

7349 ± 206 

0.25 ±0.06 

-150 ± 11 

0.82 ± 0.04 

0.14 

-0.08 ±0.03 

2.3 

0.0046 

14145 

4045 ± 100 

9255 ± 81 

0.41 ± 0.04 

—3 ± 7 

1.88 ±0.08 

0.11 

0.10 ±0.02 

5.5 

0.0541 

05691 

4640 ± 119 

6010 ±314 

0.35 ±0.06 

46 ± 12 

2.54 ±0.17 

0.32 

-0.90 ±0.14 

6.3 

0.1011 


a The columns, in order, correspond to the following quantities: (1) OGLE ID, in the usual form OGLE-BLG-RRLYR plus 
the catalogue entry number; (2) orbital period; (3) time of periastron passage, in units of BJD — 2440000; (4) eccentricity; 
(5) argument of the periastron; (6) projected semi-major axis; (7) standard deviation of the fit yj SSR/[AL — N], where 
SSR is the sum of squared residuals, N is the number of data points and M is the number of free parameters; (8) rate of 
change of the pulsation period; (9) semi-amplitude of the radial velocity K = 2tul sin ij f J or \ 3 \/1 — e 2 ; (10) mass function 
/(m) = a 3 sin 3 j/P 2 , , which is connected to the stellar masses through /(m) = m 3 sin 3 ±m s ) 2 , where turii is the 

mass of the RRL, and m s is the mass of the secondary. 


tering of stars with orbital periods around ~ 10—11 yr. 
Due to the long timescale, currently OGLE is the only sur¬ 
vey capable of discovering these variables, as high-precision 
(0.01 — 0.02 mag), long-baseline data are indispensable. In¬ 
deed, the increased baseline achieved by incorporating the 
OGLE-IV data greatly improved the quality of the fits, as 
we were able to disentangle the effects of the change in the 
RRL pulsation periods form the effect of binarity on the 
O — C diagrams, thus leading to more robust binary param¬ 
eters. Note that we have succeeded in recovering RRLs with 
Blazhko periods as short as ~ 0.5 yr, which is close to the 
limit of period detection set by the data sampling - thus 
indicating that binary RRLs with similarly short periods 
would have been found through our analysis. The apparent 
lack of RRL variables in short-period binary systems puts a 
tight constraint on the inclination (due to the wide orbits) 
of eclipsing RRL binary systems, making our method much 
superior for the detection of RRL binaries. 

Their faintness (15-17 mag in I) makes the follow-up of 
our binary candidates challenging. Still, two of the candi¬ 
dates (06498 and 09879) have higher minimum masses (cal¬ 
culated from the parameters in Table [2] for reasonable mass 
assumptions for the RRL) than the RRL itself. The secon¬ 
daries in some of these systems may be stars which have 
evolved off from the main sequence, therefore bright enough 
to contribute significantly to the light of the systems. If such 
a favorable case could be found, the mass ratio between the 
two stars could be measured through the ratio of their or¬ 
bital radial velocities. 

At the distance of the bulge, the binary candidate with 
the longest major axis (09789, ~ 5.6 AU) would have a 
projected major axis of about 0.7 mas. Gaia astrometry 
will unfortunately be hampered by the high crowd i ng and 
the fa intness of the candidates. However, ISahlmann et al.l 
(120141) have demonstrated that 120 [ias astrometry with 
FORS2@VLT is attainable on multi-year timescales, for tar¬ 
gets with similar magnitudes as here. A long-term astromet¬ 


ric follow-up program might thus be feasible, in order to 
determine the inclination of these systems. 

In this study of a subsample of OGLE LCs we have 
demonstrated that RRL stars can be detected in long-period 
binary systems, provided that high-quality, extended photo¬ 
metric datasets are used. We plan to extend our analysis 
to the whole OGLE bulge RRL dataset in the near future. 
Through long-term, dedicated monitoring using these tech¬ 
niques, and adequate follow-up, we will finally have a chance 
of determining the masses for a significant sample of RRL 
stars, and thereby directly constrain the theories of RRL 
star pulsation and evolution. 
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